The chromatin insulator CTCF regulates HPV18 transcript splicing and differentiation-dependent late gene expression.

The ubiquitous host protein, CCCTC-binding factor (CTCF), is an essential regulator of cellular transcription and functions to maintain epigenetic boundaries, stabilise chromatin loops and regulate splicing of alternative exons. We have previously demonstrated that CTCF binds to the E2 open reading frame (ORF) of human papillomavirus (HPV) 18 and functions to repress viral oncogene expression in undifferentiated keratinocytes by co-ordinating an epigenetically repressed chromatin loop within HPV episomes. Keratinocyte differentiation disrupts CTCF-dependent chromatin looping of HPV18 episomes promoting induction of enhanced viral oncogene expression. To further characterise CTCF function in HPV transcription control we utilised direct, long-read Nanopore RNA-sequencing which provides information on the structure and abundance of full-length transcripts. Nanopore analysis of primary human keratinocytes containing HPV18 episomes before and after synchronous differentiation allowed quantification of viral transcript species, including the identification of low abundance novel transcripts. Comparison of transcripts produced in wild type HPV18 genome-containing cells to those identified in CTCF-binding deficient genome-containing cells identifies CTCF as a key regulator of differentiation-dependent late promoter activation, required for efficient E1^E4 and L1 protein expression. Furthermore, our data show that CTCF binding at the E2 ORF promotes usage of the downstream weak splice donor (SD) sites SD3165 and SD3284, to the dominant E4 splice acceptor site at nucleotide 3434. These findings demonstrate that in the HPV life cycle both early and late virus transcription programmes are facilitated by recruitment of CTCF to the E2 ORF.

Complete response to neoadjuvant chemoradiotherapy in rectal cancer is associated with RAS/AKT mutations and high tumour mutational burden.

BACKGROUND: Pathological complete response (pathCR) in rectal cancer is beneficial, as up to 75% of patients do not experience regrowth of the primary tumour, but it is poorly understood. We hypothesised that the changes seen in the pre-treatment biopsies of pathCR but not seen in residual tumour after chemoradiotherapy were the determinants of responsiveness. METHODS: Two groups of patients with either complete response (pathCR group, N = 24) or no response (poor response group, N = 24) were retrieved. Pre-treatment biopsies of cancers from these patients underwent high read depth amplicon sequencing for a targeted panel, exome sequencing, methylation profiling and immunohistochemistry for DNA repair pathway proteins. RESULTS: Twenty four patients who underwent pathCR and twenty-four who underwent poor response underwent molecular characterisation. Patients in the pathCR group had significantly higher tumour mutational burden and neoantigen load, frequent copy number alterations but fewer structural variants and enrichment for driver mutations in the PI3K/AKT/mTOR signalling pathway. There were no significant differences in tumour heterogeneity as measured by MATH score. Methylation analysis demonstrated enrichment for hypomethyation in the PI3K/AKT/mTOR signalling pathway. DISCUSSION: The phenomenon of pathCR in rectal cancer may be related to immunovisibility caused by a high tumour mutational burden phenotype. Potential therapy resistance mechanisms involve the PI3K/AKT/mTOR signalling pathway, but tumour heterogeneity does not seem to play a role in resistance.

The molecular landscape of well differentiated retroperitoneal liposarcoma.

Well differentiated liposarcoma (WD-LPS) is a relatively rare tumour, with fewer than 50 cases occurring per year in the UK. These tumours are both chemotherapy- and radiotherapy-resistant and present a significant treatment challenge requiring radical surgery. Little is known of the molecular landscape of these tumours and no current targets for molecular therapy exist. We aimed to carry out a comprehensive molecular characterisation of WD-LPS via whole genome sequencing, RNA sequencing, and methylation array analysis. A recurrent mutation within exon 1 of FOXD4L3 was observed (chr9:70,918,189A>T; c.322A>T; p.Lys108Ter). Recurrent mutations were also observed in Wnt signalling, immunity, DNA repair, and hypoxia-associated genes. Recurrent amplification of HGMA2 was observed, although this was in fact part of a general amplification of the region around this gene. Recurrent gene fusions in HGMA2, SDHA, TSPAN31, and MDM2 were also observed as well as consistent rearrangements between chromosome 6 and chromosome 12. Our study has demonstrated a recurrent mutation within FOXD4L3, which shows evidence of interaction with the PAX pathway to promote tumourigenesis. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Clinical-grade whole-genome sequencing and 3' transcriptome analysis of colorectal cancer patients.

BACKGROUND: Clinical-grade whole-genome sequencing (cWGS) has the potential to become the standard of care within the clinic because of its breadth of coverage and lack of bias towards certain regions of the genome. Colorectal cancer presents a difficult treatment paradigm, with over 40% of patients presenting at diagnosis with metastatic disease. We hypothesised that cWGS coupled with 3' transcriptome analysis would give new insights into colorectal cancer. METHODS: Patients underwent PCR-free whole-genome sequencing and alignment and variant calling using a standardised pipeline to output SNVs, indels, SVs and CNAs. Additional insights into the mutational signatures and tumour biology were gained by the use of 3' RNA-seq. RESULTS: Fifty-four patients were studied in total. Driver analysis identified the Wnt pathway gene APC as the only consistently mutated driver in colorectal cancer. Alterations in the PI3K/mTOR pathways were seen as previously observed in CRC. Multiple private CNAs, SVs and gene fusions were unique to individual tumours. Approximately 30% of patients had a tumour mutational burden of > 10 mutations/Mb of DNA, suggesting suitability for immunotherapy. CONCLUSIONS: Clinical whole-genome sequencing offers a potential avenue for the identification of private genomic variation that may confer sensitivity to targeted agents and offer patients new options for targeted therapies.

Rapid, highly accurate and cost-effective open-source simultaneous complete HLA typing and phasing of class I and II alleles using nanopore sequencing.

Accurate rapid genotyping of the genes within the HLA region presents many difficulties because of the complexity of this region. Here we present the results of our proof of concept nanopore-based long read polymerase chain reaction (PCR) solution for HLA genotyping. For 15 HLA anthropology-based samples and 13 NHS Blood and Transplant derived samples 40 ng of genomic DNA underwent long-range PCR for class I and II HLA alleles. Pooled PCR products were sequenced on the Oxford Nanopore MinIoON R9.4.1 flow cell. Sequenced reads had HLA genotype assigned with HLA-LA. Called genotypes were compared with reference derived from a combination of short-read next-generation sequencing, Sanger sequence and/or single-site polymorphism (SSP) typing. For concordance, accuracy was 100%, 98.4%, 97.5% and 95.1% for the first, second, third and fourth fields, respectively, to four field accuracy where it was available, otherwise three field in 28 samples for class I calls and 17 samples for class II calls. Phasing of maternal and paternal alleles, as well as phasing based identification of runs of homozygosity, was shown successfully. Time for assay run was 8 hours and the reconstruction of HLA typing data was 15 minutes. Assay cost was £55 ($80USD)/sample. We have developed a rapid and cost-effective long-range PCR and nanopore sequencing-based assay that can genotype the genes within HLA region to up to four field accuracy, identify runs of homozygosity in HLA, reconstruct maternal and paternal haplotypes and can be scaled from multi-sample runs to a single sample.

Targeted deep sequencing of urothelial bladder cancers and associated urinary DNA: a 23-gene panel with utility for non-invasive diagnosis and risk stratification.

OBJECTIVES: To develop a focused panel of somatic mutations (SMs) present in the majority of urothelial bladder cancers (UBCs), to investigate the diagnostic and prognostic utility of this panel, and to compare the identification of SMs in urinary cell-pellet (cp)DNA and cell-free (cf)DNA as part of the development of a non-invasive clinical assay. PATIENTS AND METHODS: A panel of SMs was validated by targeted deep-sequencing of tumour DNA from 956 patients with UBC. In addition, amplicon and capture-based targeted sequencing measured mutant allele frequencies (MAFs) of SMs in 314 urine cpDNAs and 153 urine cfDNAs. The association of SMs with grade, stage, and clinical outcomes was investigated by univariate and multivariate Cox models. Concordance between SMs detected in tumour tissue and cpDNA and cfDNA was assessed. RESULTS: The panel comprised SMs in 23 genes: TERT (promoter), FGFR3, PIK3CA, TP53, ERCC2, RHOB, ERBB2, HRAS, RXRA, ELF3, CDKN1A, KRAS, KDM6A, AKT1, FBXW7, ERBB3, SF3B1, CTNNB1, BRAF, C3orf70, CREBBP, CDKN2A, and NRAS; 93.5-98.3% of UBCs of all grades and stages harboured ≥1 SM (mean: 2.5 SMs/tumour). RAS mutations were associated with better overall survival (P = 0.04). Mutations in RXRA, RHOB and TERT (promoter) were associated with shorter time to recurrence (P < 0.05). MAFs in urinary cfDNA and cpDNA were highly correlated; using a capture-based approach, >94% of tumour SMs were detected in both cpDNA and cfDNA. CONCLUSIONS: SMs are reliably detected in urinary cpDNA and cfDNA. The technical capability to identify very low MAFs is essential to reliably detect UBC, regardless of the use of cpDNA or cfDNA. This 23-gene panel shows promise for the non-invasive diagnosis and risk stratification of UBC.

Validation of epigenetic markers to identify colitis associated cancer: Results of module 1 of the ENDCAP-C study.

BACKGROUND: Chronic inflammation caused by ulcerative colitis (UC) causes a pro-neoplastic drive in the inflamed colon, leading to a markedly greater risk of invasive malignancy compared to the general population. Despite surveillance protocols, 50% of cases proceed to cancer before neoplasia is detected. The Enhanced Neoplasia Detection and Cancer Prevention in Chronic Colitis (ENDCaP-C) trial is an observational multi-centre test accuracy study to ascertain the role of molecular markers in improving the detection of dysplasia. We aimed to validate previously identified biomarkers of neoplasia in a retrospective cohort and create predictive models for later validation in a prospective cohort. METHODS: A retrospective analysis using bisulphite pyrosequencing of an 11 marker panel (SFRP1, SFRP2, SRP4, SRP5, WIF1, TUBB6, SOX7, APC1A, APC2, MINT1, RUNX3) in samples from 35 patients with cancer, 78 with dysplasia and 343 without neoplasia undergoing surveillance for UC associated neoplasia across 6 medical centres. Predictive models for UC associated cancer/dysplasia were created in the setting of neoplastic and non-neoplastic mucosa. FINDINGS: For neoplastic mucosa a five marker panel (SFRP2, SFRP4, WIF1, APC1A, APC2) was accurate in detecting pre-cancerous and invasive neoplasia (AUC = 0.83; 95% CI: 0.79, 0.88), and dysplasia (AUC = 0.88; (0.84, 0.91). For non-neoplastic mucosa a four marker panel (APC1A, SFRP4, SFRP5, SOX7) had modest accuracy (AUC = 0.68; 95% CI: 0.62,0.73) in predicting associated bowel neoplasia through the methylation signature of distant non-neoplastic colonic mucosa. INTERPRETATION: This multiplex methylation marker panel is accurate in the detection of ulcerative colitis associated dysplasia and neoplasia and is currently being validated in a prospective clinical trial. FUNDING: The ENDCAP-C study was funded by the National Institute for Health Research Efficacy and Mechanism Evaluation (EME) Programme (11/100/29).

Disruption of CTCF-YY1-dependent looping of the human papillomavirus genome activates differentiation-induced viral oncogene transcription.

The complex life cycle of oncogenic human papillomavirus (HPV) initiates in undifferentiated basal epithelial keratinocytes where expression of the E6 and E7 oncogenes is restricted. Upon epithelial differentiation, E6/E7 transcription is increased through unknown mechanisms to drive cellular proliferation required to support virus replication. We report that the chromatin-organising CCCTC-binding factor (CTCF) promotes the formation of a chromatin loop in the HPV genome that epigenetically represses viral enhancer activity controlling E6/E7 expression. CTCF-dependent looping is dependent on the expression of the CTCF-associated Yin Yang 1 (YY1) transcription factor and polycomb repressor complex (PRC) recruitment, resulting in trimethylation of histone H3 at lysine 27. We show that viral oncogene up-regulation during cellular differentiation results from YY1 down-regulation, disruption of viral genome looping, and a loss of epigenetic repression of viral enhancer activity. Our data therefore reveal a key role for CTCF-YY1-dependent looping in the HPV life cycle and identify a regulatory mechanism that could be disrupted in HPV carcinogenesis.

Evaluation of the interaction between LRRK2 and PARK16 loci in determining risk of Parkinson's disease: analysis of a large multicenter study.

A recent study MacLeod et al. has shown that an interaction between variants at the LRRK2 and PARK16 loci influences risk of development of Parkinson's disease (PD). Our study examines the proposed interaction between LRRK2 and PARK16 variants in modifying PD risk using a large multicenter series of PD patients (7715) and controls (8261) from sites participating in the Genetic Epidemiology of Parkinson's Disease Consortium. Our data does not support a strong direct interaction between LRRK2 and PARK16 variants; however, given the role of retromer and lysosomal pathways in PD, further studies are warranted.

Large-scale assessment of polyglutamine repeat expansions in Parkinson disease.

OBJECTIVES: We aim to clarify the pathogenic role of intermediate size repeat expansions of SCA2, SCA3, SCA6, and SCA17 as risk factors for idiopathic Parkinson disease (PD). METHODS: We invited researchers from the Genetic Epidemiology of Parkinson's Disease Consortium to participate in the study. There were 12,346 cases and 8,164 controls genotyped, for a total of 4 repeats within the SCA2, SCA3, SCA6, and SCA17 genes. Fixed- and random-effects models were used to estimate the summary risk estimates for the genes. We investigated between-study heterogeneity and heterogeneity between different ethnic populations. RESULTS: We did not observe any definite pathogenic repeat expansions for SCA2, SCA3, SCA6, and SCA17 genes in patients with idiopathic PD from Caucasian and Asian populations. Furthermore, overall analysis did not reveal any significant association between intermediate repeats and PD. The effect estimates (odds ratio) ranged from 0.93 to 1.01 in the overall cohort for the SCA2, SCA3, SCA6, and SCA17 loci. CONCLUSIONS: Our study did not support a major role for definite pathogenic repeat expansions in SCA2, SCA3, SCA6, and SCA17 genes for idiopathic PD. Thus, results of this large study do not support diagnostic screening of SCA2, SCA3, SCA6, and SCA17 gene repeats in the common idiopathic form of PD. Likewise, this largest multicentered study performed to date excludes the role of intermediate repeats of these genes as a risk factor for PD.

A multi-centre clinico-genetic analysis of the VPS35 gene in Parkinson disease indicates reduced penetrance for disease-associated variants.

BACKGROUND: Two recent studies identified a mutation (p.Asp620Asn) in the vacuolar protein sorting 35 gene as a cause for an autosomal dominant form of Parkinson disease . Although additional missense variants were described, their pathogenic role yet remains inconclusive. METHODS AND RESULTS: We performed the largest multi-center study to ascertain the frequency and pathogenicity of the reported vacuolar protein sorting 35 gene variants in more than 15,000 individuals worldwide. p.Asp620Asn was detected in 5 familial and 2 sporadic PD cases and not in healthy controls, p.Leu774Met in 6 cases and 1 control, p.Gly51Ser in 3 cases and 2 controls. Overall analyses did not reveal any significant increased risk for p.Leu774Met and p.Gly51Ser in our cohort. CONCLUSIONS: Our study apart from identifying the p.Asp620Asn variant in familial cases also identified it in idiopathic Parkinson disease cases, and thus provides genetic evidence for a role of p.Asp620Asn in Parkinson disease in different populations worldwide.

Investigation of c9orf72 in 4 neurodegenerative disorders.

OBJECTIVE To estimate the allele frequency of C9orf72 (G4C2) repeats in amyotrophic lateral sclerosis (ALS), frontotemporal lobar degeneration (FTLD), Alzheimer disease (AD), and Parkinson disease (PD). DESIGN The number of repeats was estimated by a 2-step genotyping strategy. For expansion carriers, we sequenced the repeat flanking regions and obtained APOE genotypes and MAPT H1/H2 haplotypes. SETTING Hospitals specializing in neurodegenerative disorders. SUBJECTS We analyzed 520 patients with FTLD, 389 patients with ALS, 424 patients with AD, 289 patients with PD, 602 controls, 18 families, and 29 patients with PD with the LRRK2 G2019S mutation. MAIN OUTCOME MEASURE The expansion frequency. RESULTS Based on a prior cutoff (>30 repeats), the expansion was detected in 9.3% of patients with ALS, 5.2% of patients with FTLD, and 0.7% of patients with PD but not in controls or patients with AD. It was significantly associated with family history of ALS or FTLD and age at onset of FTLD. Phenotype variation (ALS vs FTLD) was not associated with MAPT, APOE, or variability in the repeat flanking regions. Two patients with PD were carriers of 39 and 32 repeats with questionable pathological significance, since the 39-repeat allele does not segregate with PD. No expansion or intermediate alleles (20-29 repeats) were found among the G2019S carriers and AD cases with TAR DNA-binding protein 43-positive inclusions. Surprisingly, the frequency of the 10-repeat allele was marginally increased in all 4 neurodegenerative diseases compared with controls, indicating the presence of an unknown risk variation in the C9orf72 locus. CONCLUSIONS The C9orf72 expansion is a common cause of ALS and FTLD, but not of AD or PD. Our study raises concern about a reliable cutoff for the pathological repeat number, which is important in the utility of genetic screening.

Large-scale replication and heterogeneity in Parkinson disease genetic loci.

OBJECTIVE: Eleven genetic loci have reached genome-wide significance in a recent meta-analysis of genome-wide association studies in Parkinson disease (PD) based on populations of Caucasian descent. The extent to which these genetic effects are consistent across different populations is unknown. METHODS: Investigators from the Genetic Epidemiology of Parkinson's Disease Consortium were invited to participate in the study. A total of 11 SNPs were genotyped in 8,750 cases and 8,955 controls. Fixed as well as random effects models were used to provide the summary risk estimates for these variants. We evaluated between-study heterogeneity and heterogeneity between populations of different ancestry. RESULTS: In the overall analysis, single nucleotide polymorphisms (SNPs) in 9 loci showed significant associations with protective per-allele odds ratios of 0.78-0.87 (LAMP3, BST1, and MAPT) and susceptibility per-allele odds ratios of 1.14-1.43 (STK39, GAK, SNCA, LRRK2, SYT11, and HIP1R). For 5 of the 9 replicated SNPs there was nominally significant between-site heterogeneity in the effect sizes (I(2) estimates ranged from 39% to 48%). Subgroup analysis by ethnicity showed significantly stronger effects for the BST1 (rs11724635) in Asian vs Caucasian populations and similar effects for SNCA, LRRK2, LAMP3, HIP1R, and STK39 in Asian and Caucasian populations, while MAPT rs2942168 and SYT11 rs34372695 were monomorphic in the Asian population, highlighting the role of population-specific heterogeneity in PD. CONCLUSION: Our study allows insight to understand the distribution of newly identified genetic factors contributing to PD and shows that large-scale evaluation in diverse populations is important to understand the role of population-specific heterogeneity.

Rates of hospitalisation for influenza, respiratory syncytial virus and human metapneumovirus among infants and young children.

To inform the development of a national influenza immunisation programme and the potential role of antiviral drugs in young children, we studied 613 children aged 71 months or younger who attended Leicester Childrens' Hospital during winter 2001-2002. During periods of respiratory syncytial virus (RSV), influenza, and human metapneumovirus activity, an estimated 12.2% (95% CI: 11.4-13.1), 7.1% (6.3-7.9), and 2.5% (2.1-2.9), respectively, of all medical cases assessed in the hospital were associated with these infections. The respective rates of hospital assessments for RSV, influenza, and human metapneumovirus (HMPV) were 1042 (95% CI: 967-1021), 394 (348-443), and 223 (189-262) per 100,000 children, and for admissions were 517 (465-574), 144 (117-175), and 126 (101-156) per 100,000. Few children with influenza had a prior risk factor. Children with influenza were admitted a median of 4 days after onset of illness and none was coded at discharge as influenza. We conclude that antivirals have little role in the hospital management of children with influenza. Our data provide health economists with information to evaluate the place of universal immunisation of young children against influenza. Hospitalisation rates decreased markedly with referral age, so vaccine would need to be given in early infancy for maximum benefit.